Field-Induced Disorder Point in Non-Collinear Ising Spin Chains 
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We perform a theoretical study of a non-coUinear Ising ferrimagnetic spin chain inspired by the 
^ ' compound Co(hfac)2NITPhOMe. The basic building block of its structure contains one Cobalt 

O , ion and one organic radical each with a spin 1/2. The exchange interaction is strongly anisotropic 

■ and the corresponding axes of anisotropy have a period three helical structure. We introduce and 
solve a model Hamiltonian for this spin chain. We show that the present compound is very close 
to a so-called disorder point at which there is a massive ground state degeneracy. We predict 
the equilibrium magnetization process and discuss the impact of the degeneracy on the dynamical 
properties by using arguments based on the Glauber dynamics. 
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INTRODUCTION 



Recent advances in coordination chemistry have led to the synthesis of low-dimensional magnets with unconventional 
magnetic properties Some of them are spin chains that exhibit complex magnetization behavior under an applied 

external field. The Ising ferrimagnetic spin chain Co(hfac)2NITPhOMe (CoPhOMe in the following) is the 

first quasi-lD compound in which slow dynamics of the magnetization has been directly observed. For this system 
^ it has been shown that the very long time scale required for establishment of equilibrium preclude the experimental 
determination of the static magnetization curve at low temperature. 

This very special magnet has two kinds of spins S=l/2 : half of them are Cobalt ions Co^+ and the other half are 
organic radicals NITPhOMe. They strictly alternate along the chain and have highly anisotropic exchange interactions 
. best described by an Ising-like Hamiltonian. The chain itself describes a spiral in real space so that the local anisotropy 

■ axis have a complex pattern : see Fig.(l). The primitive cell contains three Cobalt ions alternating with three radicals. 

, Each Cobalt is related to the other Cobalts by a 120° rotation around the c-axis and has a local anisotropy axis which 
makes an angle 9 of approximately 50° with c. The effective S=l/2 spin of the Cobalt and the S=l/2 spin of the 
, radical have different values and thus the chain behave as a unidimensional ferrimagnet due to non-compensation 
^0 ■ of the magnetic moments. Only part of the magnetization curve is known because the exchange energy scale is very 
' high in terms of accessible magnetic fields. Experiments show a metamagnetic jump at low fields, much less than the 
. value expected for full saturation 

' In this paper we investigate the equilibrium, static behavior of a model Hamiltonian for CoPhOMe that we believe 
^ , encompasses all the main physical characteristics of this magnet. We show that this model Hamiltonian can be tuned 
'"^ • through a so-called disorder point 7, ^] as a function of the direction 9 of the local axes and magnitude of an applied 
^ ', external uniform magnetic field. In the case of CoPhOMe, the angle 9 is not known precisely but our estimates show 
O ■ that it is very close to the critical value 9c « 55° expected for the disorder point (cos^c = l/'s/S, the magic angle of 
NMR). 

In a ID Ising system with nearest and next-to-nearest neighbor exchange it is known that under certain con- 
ditions 1^, the spin correlation function can change from a monotonic exponential decay to a damped oscillating 
behavior, passing through zero when the temperature has a special value called the disorder point. We find that a 

■ similar phenomenon is driven by the geometrical arrangement of anisotropic magnetic centers even in the sole presence 
of nearest-neighbor interactions. The Hamiltonian for CoPhOMe is complicated by the fact that the eigenvectors of 
the g-tensors do not coincide with the axis of the chain so that it is not obviously diagonal. We show that it is possible 
to choose the quantization axis of the spins so that all the eigenvalues can be computed exactly. All thermodynamic 
properties as well as spin correlations may then be obtained by a transfer-matrix computation. We are then able to 
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Figure 1: The helical structure of the ID magnet CoPhOMe. Each anisotropic center has a local axis of anisotropy Zi. These 
axes are all tilted by a common angle d with respect to the chain axis c. The angle between the projections of Zi and Zi+i on 
the plane perpendicular to c is equal to ^ . 



show that, in a well defined range of orientation, the correlation function can switch from a monotonic to an oscillating 
behavior by application of a magnetic field. The zero correlation function condition occurs at a given temperature- 
dependent field hd(T) (disorder field). The monotonic behavior appears when the ground state is ferrimagnetic while 
the oscillating behavior is related to an antiferromagnetic arrangement of the spins. 

The disorder point with extensive ground state entropy is reminiscent of recently observed disorder effects due to 
non coUinearity in ordered systems (ex. the spin ice residual entropy in DyTi207 iSr.lO. 11]). 

In section II we give an analytical solution to the eigenvalue problem for a ring of N spins; we will show that it 
is possible to have a non trivial plateau in the zero temperature magnetization curve at an arbitrary small field he- 
In section III we study the thermodynamics of the infinite chain obtaining the geometrical disorder condition in zero 
field. In section IV we discuss the consequence of our findings on the dynamical properties. Finally section V contains 
our conclusions. 



II. ANALYTICAL SOLUTION OF THE EIGENVALUE PROBLEM 

In the spin chain compound CoPhOMe 0, 0| the zth elementary unit cell for spins consists of a pair of spins 
one-half, with Hamiltonian : 

^intraccll ^ j _ ^^^^ B • t, - ^ifi 9am (B • Z,) S^' . (1) 

In this equation, the are spin-1/2 operators describing the isotropic centers NITPhOMe and the describe the 
Cobalt 8=1/2 anisotropic ions. For the isotropic organic radicals we have giso — 2. The Ising-like coupling 0| takes 
place along the local axis of anisotropy z^. The exchange coupling between the radical and the Cobalt is strongly 
antiferromagnetic The Cobalt ions are coupled to the external field only by the spin component along the axis of 
anisotropy : the unit vector defines this direction. The restriction to a t^sl Ising coupling, the form of the g tensors 
and the fact that the anisotropic centers do not interact directly allow an analytical treatment of the problem, even 
when the local directions of anisotropy z^ are not coUinear. Our calculation is performed for the particular geometry 
in which the magnetic centers are arranged in a helical structure with a threefold periodicity (see figures [Tl^ . In the 
helical structure of CoPhOMe, the elementary unit cell of the crystal consists of three cells of the type described by 
Equation J^l, with three vectors Zi, Z2, Z3. The angle of these vectors with the chain axis c is 6* and the angle between 
their projections in the plane perpendicular to c is 27r/3 (their vector sum is along c). 

Since the value of the angle 9 is not known precisely from experiments, we will solve the problem for arbitrary 
values of 9 and obtain constraints in the case of CoPhOMe. From the study of a monomeric complex of formula 
Co(hfac)2(NITPhOMe)2 where the Cobalt ion coordination is practically the same as in CoPhOMe, it is known that 
9 is close to 50°. 

The intracell Hamiltonian Eq.Q in zero field is trivially diagonal on the basis \miso) ® \mani) if we take the 
quantization axis to be the local axis of the anisotropic center. It is important to note that this axis is a function of 
the site and is not fixed over the whole system. The t,; spin will interact also with the other nearest neighbor s^+i 
with the same Ising Hamiltonian : 



/T/intcrccU jj^^i+l '^i+l 

''•1,1+1 ~ '^''i ■ 



(2) 



Zi+2 

Figure 2: The helical structure of CoPhOMe as seen from the chain axis of symmetry. The anisotropic Cobalt ions are depicted 
as ellipsoids and the radical as spheres. 



For the sake of clarity we have labelled in (j^J the two z axes with the index of their respective cell. If we want to 
express all the operators in the same frame, we can simply choose the crystal frame. Thus Zcrys is now the chain axis 
and a perpendicular direction Xcrys can be chosen along the projection of one of the vectors, say i\. The simplest 
way to pass from the local frame to the crystal frame is to rotate the tensors ; 

/ooo\ /ooo\ 

J = , 3= 0, (3) 

V J / \ gam I 

as ordinary 3x3 matrices. If 7?.(iy9, 0, ip) is the rotation matrix parameterized by Euler angles that transforms the 
cartesian components of a general vector from the crystal frame to the frame of the ith. anisotropic center, the relative 
J* (g^) tensor is transformed according to the relation J_ Ti-if, = 0) J_ li^if, 9,1^ = 0) (a diagonal matrix in 
which the first two terms are equal is insensitive to the ip rotation). More explicitly, the transformation corresponds 
to a rotation of if about the axis of the helix and a rotation of 6 about the new y axis. In our particular geometry, 
three successive local z axes are associated to the values f — 0, |7r, |7r (see figures ITl^ . in such a way that the 
effective magnetic cell of the chain contains six spins. The Hamiltonian expressed in the crystal frame is then : 

N/2 N/2 

H^Yl U-i -J'-Si+U-X-s^ - MB X! • t, + B • • Si (4) 

where the rotated tensors J' and g"^ are periodic with period three with respect to the site index i. We take periodic 
boundary conditions : site + 1 is identified to site 1. The Lande factor of the isotropic centers is insensitive to the 
rotation. 

A. Local basis rotation 

To diagonalize the generalized Ising-like Hamiltonian, we first use the local axis as quantization axis for the Cobalt 
spin Si. Hence we can replace the spin operator sf by its eigenvalue cri/2, (Ji = ±1 and we set h = 1 everywhere. 
This is not enough to convert the Hamiltonian to a diagonal form because the other spins appear through their 
projections on two special directions : Zi as well as z.^+i in Eq.Q. However, once we have used this basis we have to 
solve a set of decoupled two-dimensional problems since the Hamiltonian is now given by : 

T^iWi}) = X! ■^['^^ + '^'+1 ^i'^'] ~ 5iso B • ti - gam O^i (B • Zi) (5) 

i 

It is an explicit function of the now classical spin configuration {<7i}. 

We see now that there is an effective magnetic field acting upon the spins of the radicals but no direct interactions 
between them : 

= I] t, • hf - i^iS 9am fT, (B • Z,), (6) 
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where the effective field is a function of the applied magnetic field as well as of the nearest Cobalt spins configuration : 



h'^^{ai,ai+i,B) = ^{ai + ai+i z^+i) - Qiso B. (7) 

This problem can then be solved immediately : the eigenenergies of each radical spin are Tih°^ /2 where Ti — ±1 and 
is the modulus of the effective field in Eq.(|7I). All eigenenergies of the spin problem are thus given by : 

E{{cr^}, {Ti}) = i ^{Tj h^^{a^,a^+l,B) - fiB gam (T^ (B • Z,)}. (8) 

i 

One has now to deal with a classical spin problem. It is straightforward to obtain the equilibrium properties of this 
model for an arbitrary orientation of the magnetic field. We will concentrate on the case where the applied field is 
along the chain axis since it reveals all the physical features of this problem. From Eq.® we see that the ground 
state spin configuration is obtained when all variables Ti are equal to -1 for any applied field. The effective field takes 
only three distinct values : 



hf(a, = +l,a,+i=+l,B) = A++ = J^/i(l + cos7)+5L(j)'+25..ojCos0, (9) 



hf(a,--l,a,+i = -l,B) = A__ = J^/i(l + cos7)+.gL(^)2-2.9,,„^cos0, (10) 



hf(a, = ±l,fT,+i=Tl,B) ^ A+_- J^i(l-cos7)+.9L(j)2. (11) 

In these equations, we have defined h = l/Ltsji?, 9 is the angle between the z-axis and the vectors z^ and 7 is the angle 
between two consecutive z^ vectors. In the umbrella structure of CoPhOMe, one has cos 7 — | cos^ ^ ~ \- To get the 
ground state we can equivalently consider the simplified energy : 

E{{(Ti}) = ^^{-hf (a„a,+i,B)-AiB<?ama,(B.zO}- (12) 

i 

With only three values for the effective field it can be written as : 

E[{<Ji\) = - Jeff OiOi+x - Bas 'Ycri + Cst, (13) 

i i 

where : 

Jcff = i{A++ + A — -2A+_} and Boff i/is 5a™ B cosS + i{A++ - A — }. (14) 

From these equation it is now clear that the Cobalt spins will adopt a ferromagnetic (F) or antiferromagnetic (AF) 
configuration in the ground state according to the value of the parameter JcS- If the ground state is AF in zero 
field, increasing the applied field will ultimately turn the system to a F state. Comparing energies, we find that this 
happens for the following condition : 

A+- - A — = |/i_B I ffam B cos 9. (15) 

(We use conventions for which hb < and the F state at large field has all ai = —1). 

In zero field, the ground state is ferromagnetic for 9 < 9^ where 9^ the value of 9 for which cos 7 = 0: cos 6*^ = 1/ V3 
hence 9c ~ 55°. It is AF for larger values of 9. Geometrically at the critical angle 9c, the units vectors ii are all 
perpendicular and cos 7 = : this is the configuration called the "magic angle" in NMR. If we have a zero- field AF 
ground state and apply a field then at some point it will be replaced by the F state : this will happen with a jump of 
the magnetization since the levels do simply cross. The location of the critical field as a function of the applied field 
and the angle 9 is displayed in Fig.|(2Jl 

N 

For 9 = 9c we have a 2 "2" degenerate ground state for all configurations of the Cobalt spins. 
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Figure 3: The critical field between F and AF ground states as a function of the angle 9 for gi^o = 2 and Qani = 9-5, 9, 8.5, 8, 7.5, 7 
from bottom to top. 




0.2 0.4 0.6 0.8 1 

h/J 

Figure 4: Ground state magnetization curve for 9 = 50°, 55°, 60°, 65°. We use Qani = 9 and giso = 2 as appropriate for 
CoPhOMe. 



B. The T = Magnetization Curve 



We now discuss the magnetization curve at zero temperature of the previous model. Prom the ingredients of the 
previous section, it is straightforward to derive simple expressions for the magnetization. In the ground state we have 
the sum of the contribution of the Cobalt spins that are always completely polarized along their local anisotropy axis 
and the contribution from the radical spins that always follow exactly the local field : 

1 1 h?^ 

{M) = Y^-HB Qani CTiZi- -Hb giso -^^^y (16) 

This expression is valid for an arbitrary spin configuration. In the F state the magnetization is nonzero only along 
the z axis and equal to : 

2 1 ^,1 Jcose + fiBQisoB 

■J^iM'') = --HBQani COS0 + -IJ.B giso ^ , (17) 
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N is the total number of spins Cobalt + radical. In the AF state, the magnetization from Eq.lQBJ is zero in zero 
applied field and again only along z in nonzero field. We find then : 



(18) 



Sample curves are displayed in Fig.Q. For F ground states the magnetization process is smooth and converges to 
a saturation value {gani cosO + f/iso)|/^s|/2 which is reached however only at i? — > oo. In the AF case, for 9 large 
enough there is a jump of the magnetization when the spins reorients themselves into the ferromagnetic state. 

From the shape of the known magnetization curve of CoPhOMe, we find that we need an angle 9 which is close to 
the critical value. From the saturation value we get an estimate of the ^-factor of the Cobalt ions once we assume 
that the g of the radical is extremely close to 2. Considering present experimental data the preferred values are 
gani ~ 9 and 6* ~ 54°. 

III. FINITE TEMPERATURE PROPERTIES 

The formula ||SJ) we have found for the eigenvalues of the closed ring can be used to obtain the finite temperature 
properties, in the thermodynamic limit, with the transfer matrix method. The analytical treatment will be given only 
for the field applied along the chain axis; in this case two successive anisotropic spins experience the same field and 
the magnetic unit cell consists of one isotropic-anisotropic pair. 



A. The a-a Transfer Matrix 



If we explicitly sum over the r coordinates we can define the a-cr transfer matrix : 

T„^cr'=exp -^l3 gam h cos 9 {a + a' )^ ^ exp ^Th°^{a,a' ,B) 

T = ±l 

and obtain the well known expression for the partition function of the problem : 



Z =Tr 
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N_ N_ 

A? +A^ 



(19) 



(20) 



where A_|_ and A_ are the eigenvalues of the matrix To-.o-' • In the thermodynamic limit the largest eigenvalue dominates 

N_ 

and Z — A_^ . The spin correlation function for — s- cx) is given by : 

K_\\k~l\ 



where the \<jy^) are the two eigenvectors of the transfer matrix H19|l. 



(21) 



B. The Zero-Field Spin Correlation Function 

The first term in H21() is simply the square expectation value of each anisotropic spin and it is different from zero 
only in non zero field. In zero field the solution to the eigenvalues problem gives : 



A± = 2[cosh(i/||) ± cosh(i/||)] and 



where : 



1 



1 



1 



2' 2' 
and in zero field v-^-^ = The a-a correlation function is : 



{CTkCri) = 



cosh(i/t|) - cosh(i/||) 



cosh(i^ll) + cosh(i/||) 



\k-l\ 



(22) 



(23) 



(24) 
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Figure 5: Dependence upon the angle 9 of the Cobalt correlation length ^ in zero field for temperatures given by J/T = 
2, 1.5, 1, 0.5 from top to bottom. 



that decays exponentially with increasing spin separation for every non zero temperature. 

For > vt^i [9 < 9c) the decay of 124|) is monotonic whereas for i/ff < v-\^i (9 > 9c) it is oscillatory with sign 
(— For 9 = 9c the expression H24(l vanishes at any temperature; this is related to the fact that the energies (jHJ 
do not depend on the a_ coordinates. The thermal fluctuations lead to the same behavior in the Ising ferromagnet 
with an antiferromagnetic next-to-nearest neighbor interaction, but not strong enough to destroy the ferromagnetic 
structure of the ground state 0, ■ Here it is the 9 angle which is the tuning parameter instead of the temperature 
that makes the system pass through a disorder point a.t 9 = 9c- From the Ea. H24|) we can also obtain an expression 
for the zero-field correlation length : 

cosh(;^ll) - cosh(;y|j^) 
cosh(j/||) -I- cosh(z^f j^) 

The factor one half comes from the fact that the correlation function Ea. (|24|l corresponds to an actual separation of 
2|fc — Z| magnetic centers. 

The Ea. (|25|l gives back the Ising model correlation length when 9 = 0. From figure lO we observe that at any 
temperature while increasing 9 the correlation length decreases until it reaches 9c- When 9 > 9c it increases again, 
but in this case the effective coupling between the spins is antiferromagnetic. For 9 = 9c it vanishes, as well as the 
correlation function, at all temperatures. 



■In 



C. Field dependence of the critical point 



It is straightforward to wonder if it is possible to have a vanishing correlation length for 9 ^ 9c using the external 
field as a tunable parameter. This happens when the smallest eigenvalue of the transfer matrix is zero (see equation 
()21|l ). The eigenvalues in presence of the field are : 



A-t = X cosh(z/||) -I- X ^ cosh(i^^j^) ± \J [Xcosh(zy||) — X ^ 00811(1^^^^)]^ + 4 cosh^(i/||). 



where : 



X = exp{-!^gamhcoat 



(26) 



(27) 



The A_ = condition is fulfilled when : 

cosh(i^ll) cosh(i/x|) = cosh^(i^| J. (28) 
This equation can be solved numerically for any 9 in order to obtain what we call the disorder field h = hd for which : 



(o-kcri) - {(Jk)^ = exp 



2\k-l\ 



= 0. 



(29) 
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Figure 6: Disorder field dependence on T for 6* = 44°, 46° , 48°, 50°, 52°, 54° from top to bottom. 
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Figure 7: Correlation length as a function of h/J = \^b\B / J for ^ = 0.5 and 6 = 52°, 53°, 54°, 6c- 

A numerical analysis shows that the Ea. H28|l has a solution only for < < 6'c. In the Ising limit {0 = 0) the equation 
(P5|l reduces to : 

cosh2(2/3J) - tanh^ {f3gisoh/2) sinh^(2/3J) = 1, (30) 

that has no solution for any h < oo. On the other hand we know that for 6 = 6c the disorder field is hd = 0. Moreover 
the Ea. H28|l is symmetric in the exchange h <-> —h. Since there is no physical reason to think the dependence on 
6 to be discontinuous this means that the disorder field already takes all the allowed values in the range < 9 < 9c, 
leaving no solution to the Eg. 128(1 for 9 > 9c- This is in agreement with the numerical results. 

In figure © the disorder field dependence on the temperature is shown for different values of 9- This dependence 
is less strong for 9 approaching 9c- For 9 = 9c, we find a line lying on the zero-field axis. 

In figure we note that the critical point appears as cusp in the correlation length vs field plot. If we consider 
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first the case of 9 < 9c, when the field is switched on, the correlation length decreases because the fluctuations of the 
spins are progressively damped by the effect of the field; in this region the nearest Cobalt spins are ferromagnetically 
correlated {{akcrk+i) ~ (^fc)^ > 0). Then ^ abruptly drops to zero for h = hd, at the critical point. For greater fields 
the fiuctuations are antiferromagnetically correlated {{uk'^k+ii — (cfc)^ < 0), even if the field orders the Cobalt ions 
in a relative ferromagnetic arrangement. It is because of this increasing order that the correlation length, a part from 
the narrow region of the cusp, smoothly goes to zero for h oo. 

In the case oi 9 = 9c the cusp lies at the origin and so the a-a fluctuations in presence of the field are always 
antiferromagnetically correlated. From this point of view, in this particular geometry, the behavior is similar to the 
one we have for 9 > 9c. It is worth noticing that for 9 > 9c the sign of the fiuctuations is always coherent with the 
structure of the zero-field ground state, while for 9 = 9c the zero-field expectation value of each spin vanishes. 



IV. STOCHASTIC DYNAMICS 



In the last section we have seen that many different behaviors of the correlation length appear by varying the 
geometrical {9) and the external parameters {T and h) of the system. The most peculiar feature of the real compound 
CoPhOMe is the slow relaxation of the magnetization on a macroscopic scale Q . This behavior has been explained 
up to now only through a modified Glauber model that was based only on a simplified Hamiltonian that did not 
contain all the complexity of Eq.l^J. From a theoretical point of view the correlation length and the slowest time 
scale are usually related by mean of a dynamical critical exponent. Thus we can easily imagine this system to show 
a very rich scenario in the long time scale dynamics too. In this section we investigate the implications of detailed 
balance on the transition functions W{ak) and W{Tk) involved in a Glauber-like approach 0. The Eq.® m zero 
field can be written as : 



N/2 
i=l 



h O-jCTi+l 



(31) 



The ratio of the two equilibrium probabilities corresponding to the configurations (ci , ri , . . . Ufc , , Uk+i , ■ ■ ■ c « , tn_ ) 
and (cri,ri, . . . (Tfc, -t^,., cTfc+i, . . ■cfn,tn_) is : 



Peq{-Tk) exp(Tfe^^l4^) l+TfeCTfeCTfe+ltanh(^^ll/lJ-) 



Peq{Tk) exp(-rfc^^n+!ln] i _ ^,.^,.^,.+1 tanh(^ 



(32) 



which must be equal, for the detailed balance condition to be satisfied, to the ratio The (|32|) tells us that even 

if the energy gap between the magnetic and the non magnetic state can be arbitrarily small for 9 ^ 9c [v]^ ~ '^Ti)' 
the energy barrier appearing in the transition probability W{Tk) will always be Ai? ~ + ^^Ti); that does not 

vanish for any 9. 

For the sake of completeness, we write the equivalent of the (|32|l for the Cobalt spins : 



Peqi-CTk) _ 1 + ^Crfe(TfcCTfe+i + Tk-l(Tk-l) tanh(i^|| - Vi^i) 
Peq{(Tk) 1 - 5Crfe(TfeCrfe+l + Tfe_ iCTfe_ i ) tanh(l^|| - i^||) 



(33) 



In this case we see that the ratio '^□''^ '^'l' — > 1 for 6* ^ 6'c. 

Even if we do not have so far a rigorous treatment for the dynamics of this system, it seems reasonable to relate 
the quantity A_B to the activation energy of the relaxation process. For our purpose it is enough to observe that 
the presence of a jump in the T = magnetization curve, occurring at relatively small fields, and the presence of a 
tunable disorder point is consistent with the observation of the slow decay of the macroscopic magnetization. 



V. CONCLUSION 



We have introduced a magnetic Hamiltonian suited to the description of the Ising-like spin chain CoPhOMe. This 
model contains two kinds of spins with very different g-factors. We have shown that by a careful choice of quantization 
axis all eigenvalues can be determined analytically if the applied field is along the axis of the chain. As a consequence 
we have obtained the magnetization curve. This curve displays a metamagnetic jump at some critical value of the 
external applied field. The position of the jump is closely related to the presence of a disorder point of the magnetic 
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model. This disorder point can be reached theoretically by tuning the angle of the Cobalt anisotropy axis with the 
chain axis. All experimental data indicate that the real chain CoPhOMe is very close to this critical value. Right 
at this disorder point there is extensive ground state degeneracy. By using ideas taken from the Glauber dynamics 
we show that this is likely to be related to the very long relaxation times observed when performing magnetization 
measurements on this chain. There is still the problem of developing a realistic model of the dynamical properties of 
this spin chain. 
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